Stretched exponential relaxation in the mode-coupling theory for the 

Kardar-Parisi-Zhang equation 
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We study the mode-coupling theory for the Kardar-Parisi-Zhang equation in the strong-coupling 
regime, focusing on the long time properties. By a saddle point analysis of the mode-coupling 
equations, we derive exact results for the correlation function in the long time limit - a limit which 
is hard to study using simulations. The correlation function at wavevector k in dimension d is found 



to behave asymptotically at time t as C(k, t) : 
A a determined constant and B a scale factor. 

PACS numbers: 05.40.-a, 64.60.Ht, 05.70Ln, 6< l 



fc d + 4 



i.35.Fx 



with 7 = (d - l)/2, 



The Kardar-Parisi-Zhang (KPZ) jl] equation is a sim- 
ple non-linear Langevin equation, proposed in 1986 as a 
coarse grained description of a growing interface. Prob- 
ably due to the fact that is the simplest generalization 
of the diffusion equation which includes a relevant non- 
linear term, the KPZ equation also arises in connection 
with many other important physical problems (the Burg- 
ers equation for one-dimensional turbulence Q] , directed 
polymers in a random medium |^-|5| etc.). 

The KPZ equation for a growing interface, described 
by a single valued height function h(x, t) on a d- 
dimensional substrate x£ SR d is: 

d t h(x,t) = vV 2 h+^(Vh) 2 + ri(xL,t). (1) 

The first term represents the surface tension forces which 
tend to smooth the interface, the second describes the 
nonlinear growth locally normal to the surface, and the 
last is a noise which mimics the stochastic nature of the 
growth process || . We choose the noise to be Gaussian, 
with zero mean and second moment (t](x, t)r](x',t')) = 
2DS d (x — x')S(t — t'). The steady state interface pro- 
file is usually described in terms of the roughness: w — 
(h 2 (x,t)) — (h(x,t)) 2 which for a system of size L be- 
haves like L x f(t/L z ), where f(x) — > const as x — ► oo 



and f(x) 



as x — ► 0, so that w grows with time 



like t x ^ z until it saturate to L x when t ~ L z . x an( i z 
are the roughness and dynamic exponent respectively. 

The phenomenology of the KPZ is well known: above 
two dimensions, there are two distinct regimes, separated 
by a critical value A c of the non-linearity coefficient. In 
the weak coupling regime (A < A c ) the non-linear term 
is irrelevant and the behavior is governed by the Gaus- 
sian (A = 0) fixed point. The KPZ in this regime is 
equivalent to the linear Edward Wilkinson equation, for 
which the exponents are known exactly x — (2 — d)/2 
and z — 2. The more challenging strong coupling regime 
(A > A c ), where the non-linearity is relevant is charac- 
terised by anomalous exponents, which are not known 
exactly in general dimension d. From the Galilean in- 



variance j^j (invariance of Eq. (Q) under an infinites- 
imal tilting of the surface) one can derive the relation 
X + z = 2, which leaves just one independent expo- 
nent. For the special case when d = 1, the existence of 
a fluctuation-dissipation theorem gives the exact results 
X = 1/2, z = 3/2. 

We take a self-consistent approach, the mode-coupling 
(MC) approximation in which in the diagrammatic 
expansion for the correlation and response function only 
diagrams which do not renormalise the three point ver- 
tex A are retained. The MC approximation has been 
remarkably successful in other areas of condensed matter 
physics, for example in the study of structural glass tran- 
sitions [Q, binary mixtures || and critical dynamics of 
magnets |^,^|. KPZ mode-coupling reproduces the exact 
values of the exponents in d = 1. Furthermore MC equa- 
tions have been shown to arise from the large iV-limit of 
an generalised iV-component KPZ equation jnj, which 
allows in principle a systematic approach to the theory 
beyond mode coupling, by expanding in 1/N. 

Our study focuses on the large time properties of the 
KPZ equation: we predict a stretched exponential decay 
for the correlation function at large times. A phenomeno- 
logical stretched exponential law was used as long ago as 
1854 to fit electronic relaxation data for a Leiden jar ca- 
pacitor jnj, and then been rediscovered many times: fit- 
ting functions involving stretched exponentials are nowa- 
days widely used in phenomenogical analysis of relax- 
ation data (examples are dielectric relaxation Jl2| and 
glassy relaxation |Q,|l3|). However, only a few analytical 
arguments fl4]] are able to reproduce stretched exponen- 
tial relaxation in complex systems. Our prediction is in 
principle amenable of numerical verification, even though 
usual numerical techniques, mainly based on simulations 
]l5[ of discrete microscopical models which belong to the 
KPZ universality class ]l6| ], are hard in the long time 
asymptotic region. 

Mode-coupling equations are coupled equation for the 
correlations and response function. The correlation and 
response function are defined in Fourier space by 
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C(k,w) = (h(k,aj)h*(k,uj)), 
G(k, W ) = <5- d (k - k>~V - ^)(^4y), 

where (•) indicate an average over r\. In the mode-coupling approximation, the correlation and response functions are 
the solutions of two coupled equations, 



G~\k, u) = Go" 1 (k, u) + A 2 / f / ^ [q • (k - q)] [q • k] G(k q, u fi)C(q, fi) , (2) 
C(k, w) = C (k, w) + £ | G(k, W ) | 2 / ** / ^ [q . (k _ q) ] 2 C (k - q, w - Q)C(q, Q) , (3) 

where Go(k, ui) = iyk 2 — ioj)" 1 is the bare response function, and Go(k, a;) = 2D | G(k, ui) | 2 . In the strong coupling 
limit, G(k, u>) and G(k, uS) take the following scaling forms 

G(k,w) = k~ z g(uj/k z ) , 
C(k,w) = kr^ +d+zS >n{uj/k z ) , 

and Eqs. (Q) and (|J) translate into the following coupled equations for the scaling functions n(x) and g{x): 

g- 1 (x) = -ix + I 1 (x), (4) 
n(x)=\g(x)\ 2 I 2 (x), (5) 

where x = uj/k z and ii(x) and ^(x) are given by 

h (x) = P Jo d9 sin d - 2 9 / °° dq cos fl(cos 9 - q)q 2z ^r~ z dy g n(y) , 

Hx) = f /* ^sin^ 2 ej °° dq(cos9~ q f q 2 Z -s r -(.d+4- z ) joo^ dyn (^v^ n(y)j 



with P = A 2 /(2 d r(^i)7r (d+3) / 2 ), r 2 = 1 + q 2 - 2qcos9. It will be convenient to write Eqs. (0,|1) as a function of time 
t, by Fourier transforming them: 



r%(t) = h(t) (6) 



77 ->s 

(t)=h(t) (7) 
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where Ii is the Fourier transform of the real part of I\ and I2 is the Fourier transform of I2 

h{t) = 2ttP /* dtfsin^ 2 J °° dqcos9{cos9- q)q 2z ^g~k (tr z )n(tq z ) , 
= -kP 11 d9 sm d - 2 6 J™ dg(cos - g )2 g 2 Z -3 r -(d+4-2z) H (^)ft(^) . 

Assuming an exponential, or stretched exponential decay for n(i), we look for an asymptotic solution for t — ► 00 of 
Eq. (Q) in the form 

n£(t) = A(Bt)* e -l Bt l f . (8) 

We will show that such a solution exists and we will determine the exponents 7 = (d — l)/2 and a = 1. Even though 
the scaling part of the correlation and response function are coupled in the two Eqs. (||) and (0), we will see that 
the specific form of g(t) is not relevant to the large t behaviour of n(t). We will proceed as follows: we evaluate the 
right hand side of Eq. (^) in the large t limit by saddle point methods. The result of this analysis shows that any 
a < 1 allows to reproduce the exponential factor. For a < z, we can then assume that the left hand side of Eq. (]?]) 
is dominated by g(0)~ 2 n(t) (see below), thus the asymptotic form of I^it) has to be matched with g(0)~ 2 n^(t). The 
exponential factor can be matched by any a < 1, but we will see that the only value of a that allows also to match 
the power law factor Vl z is a = 1. We can then proceed with a = 1 to get 7 = and 

A = g(0)- 2 4r(4z-4) 

p 2 (d-i)/2r((d - l)/2)r(2z - 2) 2 ' 1 ' 
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Let us start with the evaluation of the left hand side of Eq. (Q). Expanding | g(x) \ 2 in even powers of x around 
x = 0, we can write (n/ \ g \ 2 )(t) — a,Qn(t) + a\d 2 n{t) / 1 dt 2 + .. where ao — g(0)~ 2 . For ri(t) of the form (||) with a < z, 
the first term of the series will dominate at large t, we can therefore assume that the left hand side of Eq. (Q) has 
the asymptotic behaviour (n/ | g | 2 )(i) ~ g(0)~ 2 n(t). Let us now turn to the asymptotic analysis of hit)- First note 
that the integral hit), is symmetric in the exchange q — ► r, 9 — > cf>, where rsin0 = qsmO (see Fig 1). This allows us 
to rewrite hit) as twice the integral restricted to the region qcos9 < 1/2: 

T 2 (t)=2TrP dq d6sm d - 2 9 dg(cos 6 - q) 2 q 2z - 3 r -( d+i - 2z) n itr z ) nitq z ) , (10) 



where cos9 — max [l/2q, 1]. In hit) the function n(t) appears in the integrand with the arguments tq z , tr z . For 
large enough t, we can always safely replace n(tr 2 ) with its asymptotic form n^(ir z ), since in the new region of 
integration r > 1/2. We have to use more care with nitq z ), since such replacement is not allowed for all q, but only 
for q > Cjt x l z , where C is a large constant. We then have hit) ~ I^it) = hit) + hit), where 

hit) = 2nPf° /tl/Z dq Jg d9 sm d - 2 9icos 9- q) 2 q 2z - 3 r-( d+4 - 2z ^itr z )nitq z ), 
hit) = 2irP /~ tl/ . dq Jg d9 sin d ~ 2 0(cos 9 - q) 2 q 2z - 3 r^ d+4 - 2z ^ (ir z ) ^(^ 2 ) . 

The contribution from J\ is negligible with respect to h as t — > oo. Let us evaluate an asymptotic expression for h 
and postpone the discussion of J\ to the end. Inserting n^(t) = A iBtyl z exp(— | Bt \ a ^ z ) in hit) gives 

hit) = 2vPA 2 iBtp \ dq d9 S m d - 2 9icos9-q) 2 q 2z - 3+ ^r^ d+4 - 2z ^e- iBt ^^ +rC "\ (11) 
Jc/t 1 /* Je 

It is immediately clear that the main contribution to this integral will come from the region where /(g, 9) = q a +r(q, 9) a 
has its minimum, i.e. from the segment 9 = 0, C /t 1 ^ < q < 1/2 (for large enough t, C jt x l z < 1/2). The value of q 
which minimises / will depend on a: for a < 1, f reaches its minimum at qo = C/t 1 ^, where fiqo,0) = C a /t a / z + \ 
(l-C/i 1 /*) | Q ~ l + iC/^/^ + OiC/t 1 / 2 ). For a > 1, the minimum is realized at q = 1/2, where /(1/2, 0) = 1/2"" 1 . 
For a = 1, all q in the region C jt x l z < q < 1/2 contribute equally, giving f(q, 0) = 1. The saddle point approximation 
in the angular integral gives: 



2 (m\^T I ' Jn P -( Bt )* (9° + l 1 -?r) I 1 _ „ |2-(d+4-2z)+ 7 n 2z-3+- ( I j fl /)d-2 -(Bt) S § g|l-g|' 



7>l/i- 2,r/M-(/J/i^ 2 ,lq< > l!l > ' !! " j I (/ |- "' M 7 

'c/t 1 / 



Performing the integral over 9 in the large i limit: 



(12) 



d—1 1/2 

• ■ ■ (-) ' iBt)^-^ J f dqe -(B^(i a +\^n | i- q \* q i>, (13) 



o 



where <j) = 2z — d — 2 + 7+ (l — a) ^rp and V' = 2z — 3 + 7 — The leading contribution for large t is then given by 

r(^±i) 



d- l\ /2 



cv 



e-W(Bt)-*^ a<l 



J 2 (t)~TrPA 2 r(^—^\ (-) (i?t)^-^ <( / e-( Bt > s a=l (14) 



(Bt 



^2"^+^ e'^^iBt)'^ a > 1 



where I = J Q dqq^il — q)^ . For a < 1, the left hand side of Eq. (Q) will be dominated by a term of the form 
gi0)~ 2 n^it), thus hit) ~ must be proportional to n^it). This can only be achieved with a = 1 and 7 = 
For a < 1, the stretched exponential factor is recovered, but trying to match the power of t leads to 7 = d -^- + ip + 1, 
which can only be satisfied by the unphysical value z = 1. For 1 < a < z, the the left hand side of Eq. (Q) cannot 
be matched by h- We then conclude that a must be equal to 1, and the asymptotic solution is then given by 
Eq. (||). The coefficient A can be determined as well by observing that 7 = d -^- gives ijj = 4> = 2-z — 3, and thus 
/ = r(2z - 2) 2 /2r(4z - 4) which leads to Eq. (§). 
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Let us now go back to the integral J±. The argument 
of n{tq z ) in Ji runs between and C z . Thus we can 
obtain an upper bound on Ji by replacing n(tq z ) by its 
maximum in this region (which is given by n(0) if nit) is 
a monotonically decreasing function). An analysis simi- 
larly to the one done for J2, can then be performed for 
Ji, leading to J\ cx (Bt)*~^r e~( Bt ) * ■ Thus the contri- 
bution to /|° horn Ji is down by a factor (Bt)~i . 

The scale parameter B is just that, and remains un- 
fixed in terms of the parameters of the KPZ equation. 
In a recent work |jl7| , we proposed a scaling ansatz for 
the correlation function as z — > 2. If such an ansatz is 
correct, in a parametrization where g(0) is finite the scale 
B^ 1 on which t varies would go to zero as z approaches 
2: B^ 1 ~ (2 — 2). Thus the asymptotic solution that we 
have found here could be read as an asymptotic solution 
for all t as z approaches 2. 

In summary, we have presented an analytical deriva- 
tion of an asymptotic solution of the mode-coupling equa- 
tions for the strong coupling regime of the KPZ equation 
in the large t limit. We predict a stretched exponential 
relaxation for the correlation function, with a power law 
prefactor. We hope that this prediction will stimulate 
numerical investigations of long-time limits generally. 

The authors acknowledge the support of EPSRC under 
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